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Abstract. We consider integrals of type . . . u1^U2\ ■ ■ ■ du, with respect to 

the Haar measure on the orthogonal group. We establish several remarkable invariance 
properties satisfied by such integrals, by using combinatorial methods. We present as 
well a general formula for such integrals, as a sum of products of factorials. 



Introduction 



The computation of polynomial integrals over the orthogonal group 0„ is a key problem 
in mathematical physics. These integrals are indeed known to appear in a wealth of 
concrete situations, coming from random matrices, lattice models, combinatorics. 

These integrals are best introduced in a "rectangular way" , as follows: 

P p q 



For some previous work on the subject, mostly of asymptotic nature, see [T], [5], [6], 
[H],[9],[10],[n],|l2],[ll],[T5],[T6]. For noncommutative versions, see |3], g], [7]. For a 
motivation for the exact computation of /(a), coming from Hadamard matrices, see |2]. 

In the one-row case the above integral takes place over the "first slice" of On, known 
to be isomorphic to the sphere 5*""^, and a standard computation gives: 



In this paper we investigate the two-row case. More precisely, we are interested in the 
exact computation of the following integrals, depending on ai,bi G 2N: 



It is convenient to consider the vectors a, b with coordinates Oj, bj respectively - in this 
paper we will use bold letters to denote vectors or matrices. We make the following 






(Sai + n- 1)!! 
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normalization, where In-i denotes the polynomial integration in the above sense, over the 
group On-l- 



/ (b) = ^"-1 (^) (^) * 



The point with this normalization is that the $ quantity on the right has a number of 
remarkable symmetry properties. We have the following result: 

Theorem A. The junction $ = has the following properties: 

(1) Flipping: ^j) = <l>(^ 

(2) Compression: $(b o) = $(b ^o"')- 

(3) Transmutation: $„(^ 2^ = (1 - l/n)<|)„+2(b)- 

In this statement the main result is the first one. Quite curiously, this simple formula 
seems to resist any kind of direct geometric approach, or conceptual understanding in 
general. In what follows we will present a heavily combinatorial proof for it. 

Our second result is an exact formula for $, as a sum of products of factorials. 

Theorem B. The values of ^ are given by 

/2a\ _ jn-iy.l ^ 4^^aM {2R)\\{2S - 2R + n - 2)\\ 

\2hJ (n-2)!!^^^^ ^ ll(2r,)!(a,-r,)!(6,-r,)! {2S + n - 1)U 

where the sum is over rj = 0, 1, . . . , min(aj, and S = T^ai + R = T^Ti. 

As an illustration, let us look at the simplest case, when only one entry of a and of h is 
nonzero. After doing some standard manipulations, we obtain the following result. 

Theorem C. The joint moments of 2 orthogonal group coordinates x,y & {uij}, chosen 
in generic position (i.e. not on the same row or column), are given by 

f (n-2)!a!!/3!!(a + /3 + n-2)!! 

Jo„ (a + n - 2)!!(/3 + n - 2)!!(a + P + n- 1)!! 

for a, /3 even, and vanish if one of a, /3 is odd. 

Let us point out the fact that this kind of technical formula is quite powerful. For 
instance the n — )■ oo behavior of the above quantity can be simply obtained by using the 
Stirling formula, and we recover in this way a well-known result of Diaconis and Shahsha- 
hani stating that with n ^ oo the variables x, y are Gaussian and independent. 

Of course, there are several other interesting formulae emerging from Theorem A and 
Theorem B. We will explore them, with full explanations, in the body of the paper. 

The present paper is purely computational. For potential applications we refer to the 
papers cited in the reference list, and to the papers cited in that papers. 

Let us mention here, however, that we have in mind two kinds of applications. With 
the notation k = Eajj, these potential applications fall into two classes, as follows: 
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(1) Case "A; fixed, n — )■ oo". This is somehow the "old" problematics, coming from 
probability, random matrices, lattice models. As already explained, our exact 
results can lead to asymptotic ones simply by using the Stirling formula. 

(2) Case "n fixed, k — )■ oo" . Here the idea, needed for instance in relation with some 
difficult problems in combinatorics [2], would be to obtain, once again via the 
Stirling formula, estimates at n fixed, for classes of matrices with Eaj^ — )■ oo. 

Now back to the general problem of computing /(a), it is our hope that the present 
results might serve as a "key input" for this purpose. Indeed, most of our statements 
seem to have some natural generalizations to the multi-row case, and the problem is of 
course to check the validity of these generalized statements. We intend to investigate the 
p-row case, with p = 3, or perhaps even bigger, in a future paper. 

The paper is organized as follows: 1-2 are preliminary sections, in 3 we present an 
alternative approach to the Weingarten formula, specially designed for the two-row case, 
and in 4-6 we state and prove the main results. The final sections, 7-8, contain a number 
of consequences and refinements in the case a G M2(N), and a few concluding remarks. 

Acknowledgements. We would like to thank J. Novak for several useful discussions. 
T.B. and B.C. were supported by the ANR grants "Galoisint" and "Granma". B.C. was 
also supported by NSERC. 

1. Group integrals 

The polynomial integrals over the orthogonal group 0.„, our basic object of study, are 
best introduced in a "rectangular form" , as follows. 

Definition 1.1. Associated to any matrix a G Mpxg(N) is the integral 



with respect to the uniform measure on the orthogonal group On- 

In this definition, and also in most of the statements to follow, we prefer to be a bit 
unprecise on the meaning of the variable n. Normally in the above situation n should be 
an integer greater than p, g; but the results in section 2 below show that the integration 
over On can be given a purely formal meaning, so that n can be any kind of variable. 

For some previous work on the subject, of asymptotic nature, see [1], [5], [6], [8], [9], 
[To] . [12], [H], [15], [16]. For a motivation for the exact computation of /(a), see [2]. 

The advantage of our rectangular formulation comes from the fact that the parameters 
p, g G N effectively measure the "complexity" of the computation. For instance in the 
one-row case (p = 1), we have the following elementary, well-known result. 
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Theorem 1.2. For any ai, . . . ,aq even we have the formula 



, . {n-l)\W-\...a 

I [ai ... ag) = 



q-- 



where m!! = (m — l)(m — 3)(m — 5) . . with the product ending at 1 or 2. 

Proof. This foUows from the well-known fact that the first slice of 0„ is isomorphic to the 
real sphere S^~^. Indeed, this gives the following formula: 



/ (ai ... ttg) = a;f . . . x^" dx 



The integral on the right can be computed by using polar coordinates and the Fubini 
theorem, and we obtain the formula in the statement. See e.g. [2]. □ 

Another well-known result, of trigonometric nature as well, is as follows. 

Theorem 1.3. At n = 2 we have the formula 

'a b\ {a + d)\\{b + c)\l 



c dj (^a + b + c + d+l)U 

where e = 1 if a,b,c,d are even, e = —1 is a,b,c,d are odd, and e = otherwise. 

Proof. When computing the integral over O2, we can restrict the integration to 5*02 = S^, 
then further restrict the integration to the first quadrant. We get: 

b\ 2 r'\ 



(^^ =e.^£ (cos t)'^+'^(sin t)'+' dt 



This gives the formula in the statement. □ 
Finally, a third elementary result about /(a) is as follows. 

Theorem 1.4. The integral I{a) vanishes unless the matrix a G Mpxg(N) is "admissible" , 
in the sense that the sum on each of its rows and columns is an even number. 

Proof. This follows by multiplying the rows or columns of m by —1, and by using the basic 
invariance properties of the Haar measure on 0„. □ 

Observe in particular that in the 2x2 case, the admissible matrices are those having 
all entries even, or all entries odd. This agrees of course with Theorem 1.3. 

2. The Weingarten formula 

Our main tool for the computation of integrals over 0„ will be a combinatorial for- 
mula, whose origins go back to Weingarten's paper [15]. In this section we make a brief 
presentation of the formula, as developed in [6], and then we present a combinatorial 
interpretation of the Weingarten matrix entries, to be heavily used in what follows. 

Given a pairing vr and a multi-index i we say that "z fits into vr" if, when putting the 
indices of i on the points of tt, each string of vr connects a pair of equal indices. 
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Theorem 2.1. We have the Weingarten formula 



L 



du= ^ 5^(i)5^(j)lVfc„(7r,cr) 



where the objects on the right are as follows: 

(1) Dk is the set of pairings of {1, ... , 2k}. 

(2) The delta symbols are 1 or 0, depending on whether indices fit or not. 

(3) The Weingarten matrix is Wkn = G^n' where Gfc„(7r,(T) = 

Proof. The idea is that integrals on the left form the orthogonal projection onto Fix{u®'^^), 
which is spanned by the vectors = Sj5^(i)ei^ (g). . -^e^^^. Since the Gram matrix of these 
vectors is < C,n,^a >= Gkn{'^,cr), we obtain the formula in the statement. See [6j. □ 

As an example, the integrals of quantities of type Ui-^jj^Ui^j^Ui^j-^Ui^j^ appear as sums of 
coefficients of the Weingarten matrix W2n, which is given by: 

-1 

1 



W, 



2n 





n 


n 


n 




n 


n 


n 






n{n - l){n + 2) 

More precisely, the various consequences at = 2 can be summarized as follows. 
Proposition 2.2. We have the following results: 

(1) I{t g)=3/Kn + 2)). 

(2) I{1 l) = l/{n{n + 2)). 

(3) I{1 l) = {n + l)/{n{n-l){n + 2)). 

Proof. These results all follow from the Weingarten formula, by using the above numeric 
values for the entries of W2n'- 

4 0\ f ^uz f ^ 3(n+l)-6 3 



OJ J — — Z_."-v'M^; n{n-l)in + 2) n(n + 2) 
2 2\ f Y-T.. / r.r.^ {n + l)-2 1 

= / M11M11M12M12 = > vr2„(7r, nn) = — — — — - = — — — - 

U Uy J ^ n[n — \)\n^2) n(n + 2) 

2 0\ /" n + 1 

u\\ux\U22U22 = w^2n(nn, nn) 



zy J n{n — lj{n + 2) 

Here nn denotes the pairing of {1, 2, 3, 4} which pairs 1 with 2, and 3 with 4. 
Observe that the first and second formulae follow in fact as well from Theorem 1.2. □ 

In general, the computation of the Weingarten matrix is a quite subtle combinatorial 
problem, and the first results here go back to [15] , [6] . A quite powerful formula, which is 
however not exactly adapted to the "symmetry searching" considerations in this paper, 
was recently obtained in [5], and was further processed and clarified in [T6] . 
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The interpretation of the Weingarten matrix that we will need here is in terms of the 
0-1-2 matrices having sum 2 on each column. We call such matrices "elementary" . 

Theorem 2.3. The Weingarten matrix entries are given by 



where a e Mk{N) is the elementary matrix obtained as follows: 

(1) Label tti, . . . , tt^ the strings ofir. 

(2) Label ai, . . . ,(7^ the strings of a. 

(3) Setttij = #{r e {l,...,2A;}|r e7ii,re aj}. 

Proof. Consider the multi-indices i,j G {1, . . . , k}"^^ given by v G and jV ^ (^r, ior any 
r e (1, . . . , k}. We have ^ttKO — ^tttt' and 5(T'(j) = ^o-o-' for any pairings tt', cr', so if we 
apply the Weingarten formula to the quantity . . . Mj2fei2fc ' ^® obtain: 



= Wk„{n, a) 

The integral on the left can be written in the form /(a), for a certain matrix a. Our 
choice oii,j shows that a is the elementary matrix in the statement, and we are done. □ 

As an illustration for the above result, consider the partitions n — and a — fnlfl. 
We have i — (112233) and j — (122133), and we obtain: 



As a first apphcation of the above Weingarten methods, let us discuss the n — > oo 
behavior of /(a). Note first that with n — )■ cxo the above matrix W2n is concentrated on 
the diagonal. This property holds in fact for any k, and has the following consequence. 

Theorem 2.4. With k — T^aa we have the estimate 



Wkn{7T, a) = /(a) 







when all aij are even, and I (a) — 0{n 



-k-l 



) if not. 
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Proof. This result, known from [8j, follows from Theorem 2.1, because with n — )■ oo the 
matrix Gkn, and hence the matrix Wkn too, is concentrated on the diagonal. See [1]. □ 

We can see from the above result, and also from the explicit formulae given in section 
1, that /(a) depends in a non-trivial way on the parity of the entries aij. 

In what follows we will basically focus on the "main case" , where all entries Ojj are even 
numbers. The general case will be discussed at the end of section 7 below. 

3. Elementary expansion 

In this section and in the next few ones we investigate the integrals of type /(a), in 
the 2-row case. It is convenient to make the following normalization, where, as usual, the 
general formality comments given in the beginning of section 1 apply. 

Definition 3.1. For a, b vectors with even entries we make the normalization 

I (l) = In-, (a) (b) $ 

where denotes the integration in the sense of Definition 1.1, over the group 0„_i. 

This new quantity $ is just a normalization of the usual integral J. More precisely, by 
using the formula in Theorem 1.2 we have the following alternative definition. 

Proposition 3.2. We have the formula: 

(s\ ^ (Ea, + n-2)!!(S&, + n-2)!! /a\ 

Proof. This follows indeed from the one-row formula in Theorem 1.2. □ 

As a first, basic example, for any one-row vector a we have $(q) = /„(a)//„_i(a), and 
according to Theorem 1.2, this gives the following formula: 

(n-1)!! (Sa. + n-2)!! 
VOy (72-2)!! (Sai + n-1)!! 

The advantage of using $ instead of / comes from a number of remarkable invariance 
properties at the general level, to be established later on. 

For A;, a; G N we let fc^ = . . . A; (a; times). With this notation, we have the following 
technical version of Theorem 2.3, to be heavily used in what follows. 

Theorem 3.3. We have the "elementary expansion" formula 

V2bJ - \\ (2r.)!(a.-r.)!(6.-n)!* ^1^^ 0^"^ 2-^J 

r\...7q I — 1 

where the sum is over = 0, 1, . . . , min(aj, hi), and A = Soj, B = Tibi, R = Sr,. 
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Proof. We use the same method as in the proof of Theorem 2.3. Let us first apply the 
Weingarten formula to the integral in the statement: 

/ ^ On 

= 5,{l^^2'^)5,{l''^^ . . . q'^n'"^ . . . q"'^)WU^, a) 

<T TT 

Now let us look at a. In order for the 5^ symbol not to vanish, a must connect between 
themselves the 2ai + 26i copies of 1, the 2a2 + 262 copies of 2, and so on, up to the 2aq + 2hq 
copies of q. So, for any i G {1, . . . , g}, let us denote by 2rj G {0, 2, ... , min(2aj, 26^)} the 
number of "type a" copies of i coupled with "type 6" copies of i. 

Our claim is that when these parameters ri,...,rq are fixed, the sum on the right 
doesn't depend on a, and provides us with a decomposition of the following type: 



^(2b)^ J]iV.(a,b)/,(a,b) 



Indeed, let us label ai,...,ak the strings of <t, and consider the multi- index j G 
{1, . . . , k}^'^ given by jV G (T^j for any r G {!,..., A;}. We have Sa'{j) = ^aa' for any 
pairing cr', so by applying once again the Weingarten formula we obtain: 

/ U^h ■ ■ ■ '^lj2AU2j2A+i ■ ■ ■ U2nA+2B du ^ ^ <^7r ( l^^S^^^ ) 5^/ {j)Wkn{T^ , a') 

= ^<5,(l2^22^)5,,,iy,„(7r,a') 

Now let us look at the integral on the left. This can be written in the form I{m), for 
a certain matrix m, the procedure being simply to group together, by using exponents, 
the identical terms in the product of Wj^ 's. Now by getting back to the definition of the 
multi- index j, we can conclude that this procedure leads to the following formula: 



/ 



Ulji . . . Uij2^U2j2A+i ■ ■ ■ '^2j2A+2B dU — 1 \ ^2R qA-R 2^-^ j 



Summing up, our claim is proved, and the quantity /r(a, b) is nothing but the integral 
in the statement. That is, we have proved the following formula, where Nr{a,h) is the 
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number of pairings a as those considered above: 

/ 



2a 
2b 



^ , ^ , , . ^ 2-^-^ 0^-« 



ri...rg 



R 



Let us compute now the coefficient Nr{sL, b). This is by definition the number of pairings 
a as above, and these pairings are obtained as follows: (1) pick 2rj elements among 2aj 
elements, (2) pick 2rj elements among 26j elements, (3) couple the "type a" 2ri elements 
to the "type b" 2rj elements, (4) couple the remaining 2aj — 2rj elements, (5) couple the 
remaining 26j — 2rj elements. Thus we have: 

^ ^2a,: 



7V,(a,b) 



n 



2r, 



2k 
2r,; 



(2r,)!(2a,-2r,)!!(26,-2r,)!! 



(2ai)!(26i)!(2ri)!(2ai - 2ri)\\{2bi - 2^)!! 



1=1 



- n 



(2n)!(2a,-2r,)!(2ri)!(26,-2r,)! 

(2a0!(2b,)! 

A (2r,)!(2a, - 2r, + l)\\{2bi - 2^ + 1)!! 

Summing up, we have proved the following formula: 

(2a.)! (26.)! 



2a 
2b 



En 



l \ {2ri) ! (2ai - 2r, + 1) ! ! (26i - 2ri + 1) ! ! ' 



^2i? qA-R 2^~^ 



It remains to convert this formula in terms of $ quantities. But this can be done by 
using Proposition 3.2. By applying it twice, we get: 



2a 
2b 



^2R. 2^~^ 0^~^ 
-^2R qA-R 2^^^ 



{2A + n 


-2)! 


!(25 + n 


-2)!! 


(n-2)!!(n 


-2)1 


!n(2a.)!! 


n(2M!! 


{2A + n 


-2)! 


!(2S + n 


-2)!!^, 


{n 


-2)! 


!(n-2)!! 





2a 



'^2R 2^^^ 0^^^ 
-^2R qA-R 2^^^ 



Thus when passing to $ quantities, the only thing that happens is that the numeric 
coefficient gets divided by J|(2ai)!! J|(26j)!!. So, this coefficient becomes: 

''ri-M = ri ,,.,,L,„ n 



(2a,)!(26,)! 



\ (2ai)!!(26i)!! fj^ (2r,)!(2ai - 2ri + l)\\{2bi - 2r, + 1)1! 

' (2a, + l)!!(2b, + l)!! 

l\ {2riy.{2ai - 2n + l)!!(2&i - 2n + 1)!! 



n 



L(2n)!(a,-n)!(6i-n)! 
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Thus we have obtained the formula in the statement, and we are done. □ 
As a first consequence, we have the following result. 
Theorem 3.4. We have the "compression formula" 

valid for any vectors with even entries a, b e and c e N^. 

Proof. It is convenient to replace a, b, c with their doubles 2a, 2b, 2c. Consider now the 
elementary expansion formula for the matrix in the statement: 

Since the numeric coefficient doesn't depend on c, and the function on the right depends 
only on C = Eq, this gives the formula in the statement. □ 

We should mention that the above formula has as well a direct geometric proof. How- 
ever, this is no longer true for the "flipping formula" in section 5 below, or for the "trans- 
mutation formula" in section 6 below. So, as a general policy, in what follows we will 
simply present combinatorial proofs for all the results. The geometric methods that we 
have so far are not powerful enough, and will be rather explained in some future paper. 

4. Triangular formula 

We explore now a problematics which is somehow opposite to the "compression princi- 
ple": what happens when "extending" the original matrix (^) with a (q) component? 

In this section we present a number of technical results in this sense, which will lead to 
a number of concrete formulae, of great use for the general purposes of this paper. The 
final answer to the "extension problem" will be given later on (Theorem 6.4 below). 

Let us begin with a basic lemma. 

Lemma 4.1. We have the "basic eoctension" formula 

Kb o)~((=-+"-^*(b)-s<-"^'K^')) 

for any a, b e (2N)^, where a*^*^ — (oi, . . . , a^-i, + 2, a^+i, . . . , a^). 

Proof. By using the trivial identity Eiifj = 1, we obtain the following formula: 
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Let US translate this formula in terms of $ quantities. According to Proposition 3.2, 
we have: 

/a\ _ (Eg, + n-2)!!(S&, + n-2)!! /a 
[hj - (n-2)!!(n-2)!!na.!!n&." 



b ; (n-2)!!(n-2)!!na.!!n&.!!(a, + l) V b 
fa 2\ _ (Eg, + n)!!(S&, +n-2)!! /a 2 



Thus our above formula translates as follows: 



b 



This gives the formula in the statement. □ 
It is convenient to record as well a "recursive" version of the above result. 
Lemma 4.2. We have the "recursive extension" formula 

valid for any two vectors a, b e (2N)^, and any c e 2N. 
Proof. We use the compression formula. This gives: 

c + 2\^^fa c 2 



h J \^b 

Now if we denote the quantity on the left by K, and we apply to the quantity on the 
right the basic extension formula, we obtain: 



K 



This gives the formula of K in the statement. □ 

As a first consequence of our results, we can establish now a number of concrete for- 
mulae. The first such formula computes all the joint moments of Uu, Ui2, U2i- 

Theorem 4.3. We have the "triangular formula" 

fa c\ _ (n-1)!! {a + c + n - 2)\\{b + c + n - 2)\\ 
\b o) ~ (n-2)!! ' {c + n - 2)\\{a + b + c + n - 

valid for any a,b,c & 2N. 
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Proof. We prove this by induction over c e 2N. At c = this follows from the 1-row 
formula, so assume that this is true at c. By using Lemma 4.2, we get: 



a c + 2 
b 



1 



[a + c + n- 1)$ 



a c 
b 



(a + 1)$ 



a + 2 c 
b 



n + c — 1 

Let us call L — R the above expression. According to the recurrence, we have: 

(n-1)!! (a + c + n)!!(6 + c + n-2)!! 



R 



(a + 1 



(n-2)!! {c + ny.\{a + b + c + n-iy.\ 
n-1)!! (a + c + n)!!(6 + c + n-2)!! 



■(n-2)!! {c + n)\\{a + b + c + n + l)\\ 
Thus we obtain the following formula: 



a c + 2 
6 



(n — 


1)!! 


(n — 


2)!! 


(n — 


1)!! 


(n — 


2)!! 


(n — 


1)!! 


{n — 


2)!! 



{a + c + n)\\{b + c + n - 2) 



{c + n)ll{a + b + c + n + 1) 
{a + c + n)\\{b + c + n-2) 



{c + n)\\{a + b + c + n + 1) 
{a + c + n)\\{b + c + n)\\ 



{{a + b + c + n) - (a + 1)) 
{b + c + n-1) 



(c + n)\\{a + b + c + n + l)\\ 
Thus the formula to be proved is true at c + 2, and we are done. 



□ 



As a first observation, by combining the above formula with the compression formula 
we obtain the following result, fully generahzing Theorem 1.2. 

CorollEiry 4.4. We have the formula 



a Ci 
b 







(n-1)!! (a + Sq + n - 2)!!(6 + Eq + n - 2)!! 



(n - 2)!! (Ecj + n- 2)!!(a + 6 + + n - 1)!! 
valid for any even numbers a, b and Ci, . . . , Cg. 

Proof. This follows indeed from Theorem 4.3 and from the compression principle. 



□ 



As a second observation, at a = the triangular formula computes all the joint moments 
of ui2,U2i. To our knowledge, this quite basic result was previously unknown. 

Corollary 4.5. The joint moments of 2 orthogonal group coordinates x , y G {uij}, chosen 
in generic position (i.e. not on the same row or column), are given by 

(n-2)\a\\p\\{a + p + n-2)\\ 



I 

Jo, 



x°'y^ du 



/o„ (a + n - 2)!!(^ + n - 2)!!(q; + ^ + n - 1)!! 

for a,P even, and vanish if one of a, (3 is odd. 

Proof. By symmetry we may assume that our coordinates are x — Ui2 and y = U21, and 
the result follows from Theorem 4.3, with a — 0, c — a, b — /3. □ 
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5. The flipping principle 

In this section we state and prove the main conceptual result in this paper. This is a 
non-trivial and quite powerful, but remarkably simple symmetry result, stating that $(^) 
is invariant under the upside-down flipping of any column of (^). 

Let us begin with the case of the elementary matrices. 

Lemma 5.1. We have the formula 

(2" _ (n- 1)!! (2a + 2b + n-2)!! 
V0« 2^ " (n-2)!! ' (2a + 26 + n - 1)!! 

valid for any a, 6 e N. 

Proof. Indeed, by using the compression principle, we obtain: 

/2'^ _ /2a 0\_ /O 2a 
1^0" 27 ~ 1^ 2h) ~ \2h 

On the other hand, by applying the triangular formula, we obtain: 

/ 2a\ _ (n- 1)!! (2a + n - 2)!!(2a + 26 + n - 2)!! 
\2h y ~ (n - 2)!! ' (2a + n - 2)!!(2a + 26 + n - 1)!! 

By simplifying the fraction, we obtain the formula in the statement. □ 

Lemma 5.2. We have the "elementary flipping" formula 

/l2s 2« _ A^* 2" 0"^ 
(^l2s 0" 2^) ~ 2'^ 

valid for any s G N and any a.b.c.d E N satisfying a + b = c + d. 

Proof. We prove this result by induction over s. At s = this follows from the explicit 
formula in Lemma 5.1, because the right term there depends only on a + 6. 

So, assume that the result is true at s e N. We use the following equahty, coming from 
the triangular formula: 

/2a 2c\^ /2a 
\2b OJ \2b 2c 

Assume a >b and consider the elementary expansion of the above two quantities, where 
Kr{a, b) denotes the coefficient appearing in the elementary expansion formula: 

. /2a 2c\ A , . fl^'- 2'"+"-'' 
^\2b J " 2^ > [y^- 2^-' 

./2a 0\ ^ J. , n^A^'' 2'^"^ ^h+<^-r\ 
^\2h 2cJ =Z.^^(«'^)*(^12^ 2^+-^j 



14 



TEODOR BANICA, BENOIT COLLINS, AND JEAN-MARC SCHLENKER 



We know that the sums on the right are equal, for any a, b, c with a > b. With the 
choice b — s, this equahty becomes: 

f \<f,f O^-^'A _ ^ / ^^ 1^'' 2"^^ QS+c-r\ 

/ ^ jvy-yQ) ( ■^2r Qa+c-r 2^~'^ j ~ / j ^ryQ; S)9 I ^2r- Qa-r 2*+^"'" j 

r=0 ^ ^ r-=0 ^ ^ 



Now by the induction assumption, the first r terms of the above two sums coincide. So, 
the above equahty tells us that the last terms (r = s) of the two sums are equal: 

^25 Qa+c— s J ^ I 1^* 0" * 2*^ 

Since this equality holds for any a > s and any c, this shows that the elementary 
flipping formula holds at s, and we are done. □ 

Theorem 5.3. We have the "flipping formula" 
valid for any vectors a, b e and c, d e N^. 

Proof. Consider indeed the elementary expansion of the two quantities in the statement, 
where Kr{a., b) are the coefficients appearing in the elementary expansion formula: 

1 2b 2d) ^ ^ii^'^^'' '^ '^^^^' ^^ qa+c-r-s 2B+D-R-sj 

X nsj ij ^ ' 

^ /2a 2d\ _ v-TTz^ ^ .^z^ M ^^ A'""^'^ 2^+^-^-^ qb+c-R-sx 
\2h 2c) ~ ^ ^ ^ri[O'i,0i)^si[aj,Cj)(P [-^2r+2S qA+d-r-s 2^+'^-^-^ I 



TiSj IJ 



Our claim is that two formulae are in fact identical. Indeed, the first remark is that 
the various indices vary in the same sets. Also, since the function Kr{a, b) is symmetric 
in a, b, the numeric coefficients are the same. As for the $ terms on the left, these are 
equal as well, due to elementary flipping formula, so we are done. □ 

Observe that in the case p = q, the fiipping principle shows in fact that the function 
$(b d) is symmetric in its entries a, b, c, d e N^. We will further develop this point of 
view in section 7 below, in the "numeric" case, p — 1. 

As a partial conclusion to the results in this paper, the two-row integrals obey to 3 
general principles: compression, extension, and flipping. We will come back to these kind 
of questions at the end of the next section, with an extra invariance property. 
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6. The two-row formula 

In this section we state and prove the main resuh in this paper: a concrete formula, in 
terms of sums of products of factorials, for the arbitrary two-row integrals. 

We know from section 3 that these integrals are subject to an "elementary expansion" 
formula, so what is left to do is to compute the values of the elementary integrals. 

These values are given by the following technical result. 

Lemma 6.1. For any a,b,r we have: 

'l2r ob\ (n-1)!! (2r)!!(2a + 26 + 2r + n-2)!! 



(n-2)!! (2a + 26 + 4r + n-l)!! 

Proof. As a first observation, at r = the result follows from Lemma 5.1. 
Consider the elementary expansion formula, with a,b & N, a > b: 

f2a\ _ ^ 4^a!b! /l2- 2^-' 0^-"^ 

*l26j~2^(2r)!(a-r)!(6-r)! U'" na-r ob-r 



Qa-r 2" 



By using the "flipping principle" , this formula becomes: 

W "^(2r)!(a-r)!(6-r)! 0"+^-^'^ 



The point is that the quantity on the left is known, and this allows the computation of 
the integrals on the right. More precisely, let us introduce the following function: 

(-j^2r 2" 

Then the above equality translates into the following equation: 

^ f2a\ A A'a\b\ , , ^ , 

*W=S(2r)!(a-r)!(6-r)!»''"> + *-^'-' 

According to Theorem 1.2 and Proposition 3.2, the values on the left are given by: 

(2a\ _ (n-l)!!(2a + n-2)!!(2b + n-2)!! 
\2b J ~ (n - 2)!!(n - 2)!!(2a + 26 + n - 1)!! 

Now by taking 6 = 0, 1, 2, . . ., the above equations will succesivcly produce the values 
of '^rip) for r = 0, 1, 2, . . ., so we have here an algorithm for computing these values. 

On the other hand, a direct computation based on standard summation formulae shows 
that our system is solved by the values of </?r(a) given in the statement, namely: 

(n-1)!! (2r)!!(2a + 2r + n-2)!! 



^r{a) = (-1)' 



(n-2)!! (2a + 4r + n- 1)!! 
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Now by using one more time the "flipping principle" , the knowledge of the quantities 
(fir{a) fully recovers the general formula in the statement, and we are done. □ 

We are now in position of stating and proving the main result in this paper. 

Theorem 6.2. The 2-row integrals are given by the formula 

/2a\ _ ^ ^^aM (2i^)!!(2^ - + n - 2)!! 

W {n-2)\\^f-^^^^ ^ ll(2r,)!(a,-r,)!(6.-n)!' (25 + n-l)!! 

where the sum is over = 0, 1, . . . , min(aj, hi), and S = Ea^ + S6j, R = Er^. 

Proof. This follows from the elementary expansion formula, by plugging in the explicit 
values for the elementary integrals, that we found in Lemma 6.1. □ 

As a first remark, all the results in the previous sections, and in particular the com- 
pression, extension, and flipping principles, can be deduced from the above formula. 

We would like to present now another invariance principle, which is somehow of different 
nature, because it involves a "transmutation" of the n variable. The strange nature of 
this principle comes as well from the fact that is not clear how to obtain it directly. 

Lemma 6.3. We have the "basic transmutation" formula 

a 2\ n - 1 ^ /a 



valid for any two vectors a, b e (2N)^. 

Proof. It is convenient to replace a, b by their doubles 2a, 2b. Consider the elementary 
expansion of our two quantities, where K denotes as usual the numeric coefficients: 

^ /2a 2\ _ T-r ^ , , 2^-^+1 0^-«\ 

1 2b / ~ 11 1 1^^ qA-r+1 2^-^) 

^ ' ri...rq i ^ ' 

^ /2a\ V- TT . , ^^ A'"" S^""" 0^-^\ 

"+^i2bj ^ 2-^ [[J^ri[ai,0i)<Pn+2 {■^2R qA-R 2^-^) 
^ ' r\...rq i ^ ' 

In these two formulae the sums are over the same indices, and the numeric coefficients 
are the same. So, it is enough to prove that we have the following equality: 

/l^R 2^-^+1 0^~^\ n — 1 /l^^ 2^~^ 0^~^ 

n 



■^2R qA-r+1 2^~^ J ~ n +^ yi"^^ 0"^"^ 2^~^ 

That is, we just have to prove the basic transmutation formula for the elementary 
matrices. But this follows from the explicit formula in Lemma 6.1, and we are done. □ 
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Theorem 6.4. We have the "transmutation formula" 

/a c\_(n-l)!! (SQ + n-2)!! 
^^[h Oj (n-2)!!'(Eci + n-l)!! "-^''^Vb 
valid for any vectors a, b e (2N)^ and c e (2N)p. 

Proof. First, by using the compression principle, we just have to prove the above formula 
at p = 1. That is, we have to prove the following formula: 

a c\ _ (n-1)!! (c + n-2)!! /a' 
b o) ~ (n-2)!! ■ (c + n-1)!! 

But this follows from the basic transmutation formula. Indeed, by applying this formula 
c times, and by using the compression principle, we obtain: 

a 2c\ _ n-1 n + 1 n + 2c - 3 /a' 
b oj - "^■;w^---n + 2c-2 "^''^ [h, 
(n + 2c-2)!! (n-1)!! 



a 

n-2)!! (n + 2c- 1)!! 

a 



*n+2c 

- ±J11 

(n-1)!! (n + 2c-2)!! 



(n-2)!! (n + 2c- 1)!! V^y 
This gives the formula in the statement. □ 

7. The 2x2 case 

In this section we discuss the case of the 2x2 matrices a e M2 (N) , with a number of 
refinements of the results in the previous sections. As we will see right away, there is some 
"magic" in the 2x2 case, waiting to be fully discovered, and conceptually understood. 
The present section should be rather regarded as an introduction to the 2x2 problematics. 

The "magic" comes from the following result. 



a c 
b d 



Theorem 7.1. The function 

^' " {a + d + n- 2)!!(6 + c + n - 2)!! 

is symmetric in a, b,c,d. 

Proof. This follows from the flipping formula for the 2x2 matrices. Indeed, by using the 
conversion formula in Proposition 3.2, we obtain: 

(n-2)!!(n-2)!!a!!6!!c!!d!!$ ^ 
f{a,b,c,d) = + ^ + ^ _ 2)!!(& + c + n - 2)!!(a + c + n - 2)!!(fo + d + n - 2)!! 
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The flipping principle tells us that the $ quantity is symmetric in c, d. Now since the 
coefficient is symmetric as well in c, rf, we conclude that / is symmetric in c, d. Together 
with the standard fact that / is symmetric in a, d, and also in 6, c, this gives the result. □ 

The problem now is to find a formula for /(a, 6, c, (i), as a "sum of products of factorials, 
symmetric in a, 6, c, d" . Let us introduce the following basic quantities. 

Definition 7.2. Associated to a, 6, c, are the following quantities: 

(1) Sk{a,b,c,d) = {a + k){b + k){c + k){d + k). 

(2) Ffc(a, 6, c, rf) = (a + A;)!!(6 + A;)!!(c + k)\\{d + k)U. 

(3) Pk{a, b,c,d) = {a + b + c + d + k)\\. 

Observe that all the above quantities are symmetric in a, b, c, d. 

Theorem 7.3. We have the following formula: 

{n-2)\{n-2)U Fq 



f{a,b, c,0) 



n-l 



Proof. From Theorem 4.3 we get: 

'a c\ (n-l)!! (a + c + n - 2)U{b + c + n - 2) 



b Oj (n-2)!! {c + n - 2)\\{a + b + c + n - 1) 
By using the above conversion formula between / and we obtain: 

, (n-2)!a!!fe!!c!! 

f{a,b, c,0) 



{a + n- 2)!!(6 + n - 2)!!(c + n - 2)!!(a + 6 + c + n - 1)!! 
This gives the formula in the statement. □ 

In the general case, Theorem 6.2 provides of course a concrete formula for /(a, b, c, d), 
as a double sum of products of factorials. However, since that formula is not obviously 
symmetric in a, b, c, d, there is definitely room here for some improvements. 

Instead of getting into this subject, which is a bit away from the purposes of this paper, 
let us just write down the final formula in the 2x2 case, in the form of a conjecture. 

Conjecture 7.4. For a, b, c, d even we have the formula 

, (n-2)!(n-2)!! + At - ^ (n + 2r - A\ S^S.^ ■ ■ ■ S.2r+2 
f{a,b,c,d) = — > — 

where the sum is actually finite, stopping at I = min(a, b,c,d)/2. 

As already mentioned, this formula should follow from Theorem 6.2, by doing some 
summation work. While the techniques here, such as Sister Celine's method, don't lack 
(see [13]), we would rather keep this technical work for one of our future papers. 

Let us also mention that the above formula was actually obtained at an early stage of 
the present work, and comes with heavy computer evidence. 



■r-3 



ON POLYNOMIAL INTEGRALS OVER THE ORTHOGONAL GROUP 



19 



Finally, in the case where a, b, c, d are odd, the situation is quite similar. An analogue 
of Theorem 7.1 holds, in the sense that the function / given by the formula there can be 
shown to be symmetric in a, b, c, d. We have the following conjectural formula: 

Conjecture 7.5. For a, b, c, d odd we have the formula 
f( h d\- (^-^)'^" Fi ^ n + Ar-l /n + 2r - 3\ ■ ■ ■ g^2r+i 

^ ' ^ Pn-l >n-l ^ - 1) (n - 3) V 2r + 1 S„5„+2 . . . 5^+2.-2 

where the sum is actually finite, stopping at I = {min{a,b,c,d) — l)/2. 

Once again, this is a statement coming with heavy computer evidence. As regarding a 
potential proof, this should come from a suitable extension of Theorem 6.2 to the "odd" 
case. But getting into this technical subject is beyond the purposes of this paper. 

As a last remark, it is not clear how to unify the above "odd" formula with the "even" 
one from Conjecture 7.4. A more elaborated definition for the quantities S, F, P, and 
perhaps even for the double factorials themselves, is probably needed here. 

8. Concluding remarks 

We have seen in this paper that the two-row integrals of type /(a) enjoy remarkable 
symmetry properties, which can be effectively used for their exact computation. 

It is our hope that the present results will substantially contribute to the further de- 
velopment of the general study of integrals of type /(a). As explained in [2], the question 
of computing exactly these integrals, which perhaps lacks a bit of motivation in the gen- 
eral context of "probabilistic" mathematical physics, where the n — )■ 00 limit is usually 
the correct quantity to look at, is however a very interesting one in the context of pure 
mathematics, due to its relation with the Hadamard conjecture. So, our hope is that the 
present work will be a useful adding to the lineup of recent papers [2], [5], [9], [1], [16]. 

Finally, let us mention that there are some interesting "noncommutative analogues" of 
the problems investigated in the present paper. The free analogue of the hyperspherical 
law was recently found in [5], and some further occurences of this law come from [1], [7]. 
However, many questions are still open. For instance computing the "exact correlation" 
between 2 coordinates following this law remains a remarkably difficult, open problem. 
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